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ON INTEEMEDIATE ORBITS. 

By Db. G. W. Hill, West Nyack, N. Y. 

The assumption of the Keplerian ellipse as the first approximation to the 
motion of a planet leads to some inconveniences, the worst of which, perhaps, 
is that the mean longitude at the epoch suffers a perturbation proportional to 
the time, and thus the mean motion in longitude is not the same in the perturbed 
and unperturbed orbits. In the case of Saturn this perturbation amounts to 
110" per year. This inconvenience would be avoided if, developing the per- 
turbative function in a series proceeding according to the cosines of multiples 
of the angle H contained at the sun by the radii of the planets thus 

R^^f^ (r, /) + ^1 {r, Z) cos II -\- f^ {r, r') cos 211 -f . . . , 

we should suppose (p^ (r, r') annexed to the potential of the attraction exerted 
by the sun and the resulting differential equations integrated and the expres- 
sions of the coordinates thus deduced regarded as the first approximation. 
Then the following approximations could be obtained by the method of varia- 
tion of the elements or otherwise by attributing to R the value 

R ^ <Pi {r, /) cos II + fii')', r') cos 211 + . . . . 

In this way of approaching the question, the first approximation to the 
motion of a system of planets revolving about their central body would involve 
a potential function containing the radii of the planets alone, without their 
longitudes or latitudes ; thus the force acting on each planet would be directed 
along its radius, but would be a function of all the radii. The orbits are then 
all plane curves whose planes pass through the central body, and thus we know 
precisely what functions the latitudes and the reduced longitudes are of the 
orbit longitudes. Also there is equable description of areas by the radii, and 
when the latter are known functions of the time the determination of the orbit 
longitudes is rediiced to quadratures. Thus the question is narrowed to the 
finding of the radii as functions of the time. 

Let r, r', . . . v, v', . . . denote the radii and orbit longitudes of a system of 
planets, m, m', .. . their masses and JI the mass of the central body ; // = M-\- m. , 
/i' = Jl/" -|- »«',... ; R,R', . . . the perturbative functions for m, m', .... Then 
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the differential equations for determining the radii are 
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But, h, h', . . . being constants, we have 

dv h dv' h' 

di~^' 'di~¥^' ■■■■ 

The functions R, R', . . . are to be reduced to 

R = m'f (r, /) + m"f {r, /') -|- . . . , 
R' = m<p (/, r) + »'t"f (r', /') + ..., 



where we have the relation f {x, y) = f (y, x) . To evaluate the function f, let 
the symbol SK denote the operation of taking the arithmetico-geometrical mean, 
then 

W(7', /) = - = ^ 

^^ ' ^ 9M (/ -f 7-, / — r) sm(r-', x/r'^ — r=) ' 

where it is supposed that r' ^ r. The two quantities under the symbol SW may 
be brought still nearer to equality by continuing the operation, and writing 
!,/?•' = a , ],V/^ — r' r=: y9 , we have 

1 1 



(p{r,r') = 



m [i («^ + /?), «/9] SK [J (a+^/, i/^ a^ («2 + ^)] 



And if r and / are not too near each other a sufficiently approximate value is 

/ '\ 21^2 
(D (r, r) = . ' . 

If we write 

Q /JLm fi'm' 1 mh^ 1 m'h'^ 

-|- mm'ip (r, r') + mir)!'f {r, /') + m'in"(p (/, /')+..., 

the differential equations determining r, /, . . . become 
d'^r _9ii ,d^r' __9Q 

"^ de~ dr' "^ 'dF~ dr" •••• 
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From which we derive the integral equation 



m 



^. + rr."^+... = H^J+C), 



C being the arbitrary constant. The left member of this cannot be a negative 
quantity. Consequently, if we construct in a space of n dimensions, n being 
the number of planets, the surface whose equation is 

/^+ C=0, 

the representative point P, whose coordinates are the values of the variables 
»', /,..., must lie on the positive side of this surface. In such a system of 
planets as our solar system composed of the eight major planets, at least one 
fold of this surface is closed ; and it is within this that the representative point 
P always falls. Thus, in this case, we are able to set definite inferior and 
superior limits to >■,/,..., which, however, in general are not the rainhna- 
minimorum, or maxirna-maximorum values of the variables, since the point P 
only attains the surface when all the planets are together on the lines of their 
apsides. 

A particular solution of the foregoing system of differential equations can 
be obtained in the following way : Put 

^A — ^A — o 

9r~ ' 3/ " ' ■ ■ ■ ' 

and solve these equations, regarding r, r', . . . as the unknowns. Let 

r = a , r' --—a! , . . . 

be a system of values satisfying them. Then the last equations evidently sat- 
isfy the differential equations. When this process is applied to the solar 
system the values found for the a do not greatly differ from the mean distances 
of the several planets from the sun. The excessive smallness of the in rela- 
tively to M brings it about that there is but one real solution ; and quite ap- 
proximately, by neglecting in Q all the terms dependent on the functions tp, we 
get 

fi n 

By substituting these values in the terms arising from these functions we can 
arrive at more exact values for the a ; and this process can be repeated as 
often as is deemed necessary. Evidently all this comes to the adjustment of 
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the planets at such distances from the sun that for the given values of the w. 
and the h the centripetal force may be equal to the centrifugal force in each 
case, thus rendering it possible that all the orbits may be circles. 

Desiring now to find inferior and superior limits to the particular variable 
/'^, in the surface ,G -f C ^ we suppose this variable to be a function of all 
the rest, and thus that when it arrives at a maximum or minimum value we have 

5r "~ ' 5/ ~ ' ■ ■ ■ ' 

except that ^s—j^^ = 1 . In order to find limits for the values of t^^ we take the 
group of equations 

~9^~ ' 9r'~ ' -•-' 

SQ 
and in it replace the equation — ^^ =r by the equation ii -\- C = . Solving 

these equations, regarding r, /, . . . as the unknowns, we shall always, in cases 
like our solar system, be able to find two sets of corresponding real values for 
these quantities such that the two values of /" shall contain between them the 
observed value of this radius from which the value of C was derived. And no 
other roots can be found which satisfy this condition. In ordinary cases these 
roots are non-multiple and there is but one solution to the problem. 

For the sake of illustration, limiting ourselves to two planets, we may take 
a hypothetical case, suggested by Jupiter and Saturn. Assume as the values 
of the planetary masses 

1 ...., 1 



1047.355' 3601.6' 

From Leverrier's Tables for the time 1875, Jan. 1.0, we get 
log r = 0.7367630 , log / = 0.9953504 ; 
and, with the Julian year as the unit of time, 

log ~ = 8.0841634 n , log ^' = 8.8776073 . 

From the same source the values of the osculating elements necessary for the 
computations are 

log a = 0.7162505, log e = 8.6865018, log a' = 0.9794334, log e' = 8.7291462. 
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Accommodated to the units of length and time we employ, 
log// =1.5967563, log//' = 1.5964658, log wi = 8.5762479, logm' = 8.0520753. 
log h^ = 2.3119809, log h'^ = 2.5746499. 

When these numerical values are substituted in the equation Q -\- C = 0, it is 
found that C = — 0.166513246. Employing brackets to denote numbers cor- 
responding to common logarithms, the equation i? -f 6'= becomes in this 
case 

[0.1730042] r-'— [0.5871988] r-2 + [9.6485411] r'-' — [0.3256952]/-^ 



+ [6.6283232] r-' 4> 



;^] = 0.166513246. 



/j 



And the equations -r— = 0, -— = become 

^ dr ' 9r 

2 [0.5871988] r-' — [0.1730042] + [6.6283232] ^fl-,] = , 

2 [0.3256952] r'-^ — [9.6485411] — [6.6283232] \</>[^ + ^ f f^l 1 = , 

where </> denotes the same function as Laplace's bf is of «. 

The solution of the first and third of these equations gives inferior and 
superior limits for r, the solution of the first and second the same for /. Em- 
ploying the tentative process we obtain as the corresponding values of the 
logarithms of the radii in the four solutions 

log r = 0.6926819 , 0.7390339, 0.7152497, 0.7152334; 

log/ = 0.9776575, 0.9776183, 0.9238410, 1.0390527. 

The solution of the second and third equations gives the values of the 
radii for which the planets could describe uniformly circles. These are 

log a = 0.7152396 , log a' = 0.9776403 . 

If we assume that the maximum eccentricities of Jupiter and Saturn are 
respectively 0.0594902 and 0.0845677, the minimum and maximum values of 
log r and log / are severally 

log r = 0.6895907, log r = 0.7413343, log/ = 0.9411223, log / = 1.0147527. 

The differences between the former values and these are, of course, due to 
the neglect of the terms of E involving the angle If. 
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For the purpose of graphically exhibiting the numerical results just ob- 
tained we may suppose a system of rectangular axes to be drawn with r and 
r' as coordinates. The four sets of corresponding values of these variables 
give four points on the oval which constitutes one of the branches of the curve 
having Q -\- (7 = as its equation. These suffice to roughly indicate the course 
of the oval within which the representative point P, exhibiting the simulta- 
neous values of the radii, always keeps. The position for which both radii are 
constant is well towards the centre of the oval. If we draw the right line 
bisecting the right angle at the origin formed by the axes of coordinates, the 
intersection of this line or its nonintersection with the oval will show whether 
there is intrusion or not of the spheres of the planets on each other. In the 
case here treated there is no intrusion ; each planet maintains its character as 
exterior or interior ; consequently the distance between them can never vanish, 
and stability of motion is assured. 

The given example shows that there must exist a large domain in which it 
is possible to expand the function i? in a series of powers and products of the 
deviations of the radii from their mean values. Here we will limit the exposi- 
tion to the case of two planets, as the formulae given can be readily extended 
when there are more. Let us put 

r = a -\- — 7= , / = a' + 



where the second terms are small relatively to the first. The differential equa- 
tions then become 

de~ 3x ' de ~ 3x'' 

In developing Q in a series of ascending powers and products of x and x' we 
will stop with terms of three dimensions. The constant term can be omitted 
as its retention would only serve to modify the value of the arbitrary constant 

10 10 

G. As a and a' are corresponding roots of the equations ^r^ = 0, ^^, = 0, the 

terms of one dimension vanish. Also the term in axtf can be supposed absent, 
for, if present, it could be removed by a linear and orthogonal transformation 
of variables. We write therefore 

i^ = — i aV — i a'V + icx' + bx'x' + I'mx"" + ^ c'«'' , 

where a, a', h, h', c, c' are constants.* 

* The reader is asked not to confound these a and a' with the quantities previously denoted by 

the same letters. 
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Adopting two constants n and n' at present left indeterminate, but to be 
determined hereafter so as to fulfil certain conditions, the differential equations 
can be written 

-rp -\- n^x == (w^ — a^) X -\- c!>? + ^bxx' + b'x'^ , 

^ + 7«'V = (r,,'2 — a'^)x' + c'a;'^ + 'ib'xx' + M . 

Before integrating these it will be well to ascertain superior and inferior 
limits for x and x'. Those for x are determined by the equations 

C — iaV — i a V^ +icx^+ bx'x' + Vxx^ + \ cV^ = , 

— a'V + dx"^ + Wxai + Sa^ = . 

If we put )y = ± ^ , this will be a small quantity of the same order as x 

and the limiting values of x can be developed in an infinite series of powers 
of Tj. Thus 

where we take rj positively for the superior limit and negatively for the inferior. 
The corresponding value of of can be derived from the equation 

, /o77 

The limits for x', by putting ;; = + ' y , can be obtained from those of x by 

removing the accent from the symbols which there have it and applying it to 
those which are destitute of it. 

By neglecting the right-hand members of the differential equations as last 
written their integrals are simply 

x ^ e cos (f , x' = e' cos <p' , 

e and e' being arbitrary constants and f and <p' are arguments increasing pro- 

portionably to t so that -^ =^n , -^ = n' . By substituting these expressions 
eft cct 

for X and x', 

ca^ -\- '2bxx' + b'x'^ = ^ {ce^ + b'e^) -\-\c^ cos 1<p + bee' cos (tf + tp') 
+ bed co8(f — ¥'') + i ^'^'^ ^'^^ ^f' ' 
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c'a-"'^ + 25'a;a;' + So:^ = ^ {ce"^ -^-hi) -\-\ cV^ cos 2^' + h'ed cos {tp + f ') 
+ Ved cos (^ — *p') -\- \ ^^ cos 2^ . 

Whence it follows that x contains the additional terms of the second order 
with respect to e and e' 



h (-^ + *'^'^) - 6^^ '' ^^« 2^ + ^?^+ nj 



^=oZ2 io^' + *'«") - «;r2 ''^ ««» 2^ + -, 7— T-'^T^ cos {<p + ^') 



+ ^ 7 TB cos (f — f ) + o -1 1-^ cos 2y7 , 

and af terms which are obtained from those of x by simply interchanging the 
accents. 

Pushing the approximation to terms of the third order, a^ is found to con- 
tain the following terms of this order : 

x' = [« ce' + b'e'^ -2 cos <p + ^^, __ ^„ cos <f' - g— , cos 3^ 

+ -2 7 ^ rTi COS (2^ + ^C') + -2 ^ TTj COS (2^ — ^ ) 

, 1 Vee"^ , , n /\ , t b'ee'^ , „ « 

+ 2;F=:4^^°'(^ + ^^) + 2,?^^4^'^^'(^ -^^>- 

The similar terms of a;'^ are obtained from those of a? by interchanging the 
accents. The terms of the third order of ««' are 

, rl c'e'^ + b^ , b(e"-ie')-\ , Flce'+b'e'' , 5'(«'-iOl ' 

-[21^ + M 2? -^ (^ + V) + [2^, - i] £cos(, -2,'). 

The coefficients w^ — a^ and w'^ — a'^ are of the second order, hence limiting 
ourselves to terms of the third order 

(n^ — a'')x = (n' — a^) e cos <p , {n'^ — a'^) a/ = (n'^ — a'^) e' cos (p' . 
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The constants e and e' being arbitrary, the coefficients of cos ^ in a? and of 
cos (p in «' need receive no corrections in the following approximations. Hence 
the coefficient of cos ^ in the right member of the first differential equation 
should vanish, as also that of cos <f' in the second. This gives us the two 
equations 

Vc , hd , W , W l.,2 Q 



\_n^ ^ ii'^ ^ in' — n" ^ 4.n'^ — «==] 






of which the second is obtained from the first by interchange of accents. To 
the degree of approximation to which they are pushed they serve to determine 
the values of n and rt as functions of the constants a, a', b, b\ c, c', e, e'. As 
the next terms which would appear in these equations through the following 
approximations are of the fourth order with reference to e and e', the errors of 
n and n', as determined from these equations, are of the same order. 

The remaining terms of the third order in the right member of the first 
differential equation are 

tbo , I be . bb' , , bb' ~\ Ji , /o , ,> 

2^-+^ + ^ ^ + 2^,-Tl? + * il^^^J ^' '°^ (2^ + ^ ^ 

+ bJ^^ - ^ I + ^^>^' - * iJ^l^] ''' ««« (^^ - ^') 

[, b'c , , be' P b'^ ~\ ,2 f n >\ 

* 4;^7? + V^ - 2^^7^7?-^ - W'-^7?J^^ "^^^ (^ - 2^ ^ • 

The similar terms of the right member of the second differential equation are 
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obtained by interchange of accents in the preceding expression. The terms 
of the third order in x are then 

2nn' - ;r2 ~ T ^ + 2«Z^ ^ * 4^2 _ ^^'^J 3„2 _ inn' + n'' ^^^ "^^ ^^ 
i 4^-i^^ + ^ ;?^ + 2nn' + n'' + 2n«' + n^J 4nn' + in'' *'''' ^^ + ^^'^ 

r, 5'c ,1 ;'6-' p b'^ -| g«'' / _9 '\ 

L* 4;i'2 _ ^2 "^ ^ n-^ 2«w' — n'^ Inn' — w^J Ann' — 4«'2 "^ ^^ ^*P > ' 

The similar terms in x' are obtained from these by interchange of accents. 

On many occasions it will be interesting to know the equivalent of C, the 

arbitrary constant attached to the integral equation, in terms of a, a', h, h', c, c', 

da? 
e, e'. To obtain this we compute the non-periodic terms of a^, -^-^ , a^, a?x' . 

dx'^ 
The non-periodic terms of x'^, ~j^ , «'', x'^x can be derived from these by inter- 
change of accents. Limiting ourselves to these terms we have 

«?«' _ 1 n V -f ^ e* + 4»^ - dnhi" + n" ^, ,^ n" .„ 

W~^ + 18«^ + w'^ (4n^ - n'2)^ e e + * ^^^,, _ ^^,^, e , 

8^2 — 3w'2 



^ n^(4n^ — n^) L4n^ — n^ 4w^J 

Substituting these values in the integral equation and eliminating n' and w'^ by 
means of the two equations which determine them, we find 

6_i(ae+aO L 16«'^(4«^-«'7^ '^ +144^J' 

, f Sa' — 18a V^ -\- 5a'* ,2 8ft'* — 1 8aV^ + 5a^ ^.^ 3 b'c „ '^O 2-2 
+ L 2a'2(4a^-0' 2a2 (4a'2 - a'O' * «'^ ^a'^J^^ 

r 96a'* — 68rt%'- + 9a* ^ , 37 c'^-l ,, 
L 16a^(4a'2 — a'^f '^144a^J* ' 
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It is necessary to notice the effect of the vanishing of the divisors intro- 
duced by the integration. In the terms of the second order if w = 2n' the 
coefficient of cos 2f ' in x and the coefficient of cos {<p — <p') in x' apparently 
become infinite. Also, if n' = 2n, the coefficient of cos (<p — <p') in x and the 
coefficient of cos 2^ in x' are in like case. This evidently is the analytical 
warning that in these cases the arguments 2f , 2<f', f — f' cannot be distin- 
guished from f and f. Thus the coefficients of their cosines should have 
some indetermination. But let us see what conditions the quantities a, a', h, b', 
c, c', e, e' must satisfy that either of these relations may liave place. We shall 
not allow that any of the quantities a, a', n, n' can be infinite. Then, from the 
equations determining n and n\ it is plain that, in the first case, either 5' ^ 
or 2e^ — e'2 := , and in the second case, either 5 = , or e^ — 2«'^ = 0. Thus 
all the mentioned coefficients in x and x', instead of becoming infinite in these 
cases, really take the form ^. 

The particular supposition of e = 0, e' = makes the first approximation 
vanish, and it is plain from the differential equations that a; = 0, «' ^ form 
a particular solution of them in which n and n' are quite indeterminate. This 
case therefore does not demand further consideration. But suppose that, in 
the first case, J' := or, in the second, 5 = 0. Then, in Q, the term h'xx'^ or 
the term ha?x' disappears, and the terms in x and x' having indeterminate co- 
efficients should also disappear. However, the relation n = %i' or w' := 2n , 
without being exactly fulfilled, may be so very nearly. Let us take the first. 
If, in the equations determining n and n' we make n^ = 4w'^ and w'^ = a'^ except 
in the first terms and those which involve the divisor 4w'^ — w^ , we shall have 
quadratic equations for the determination of ■n? and ri^. Thus, putting 

«'^ = a'^ - (I c'^ + i h'^) ^, - (i b'c -f he' + ^\ b') ^, , 

^ = ia^~a'' + (I c'^ + ib'^- tV b'c - i be' - ^tt b') J. 
+ (ib'c + bc'-^^^\P-^\(^)^^ 
the solution of these quadratic equations will give 



n^ = 4«'2 ± 4i//32 — b'\2e^ — e'^) , 



In"" - »,2 = — 2j8 + 2i/''j8^ — b"" {2e^ — e'^) . 

The ambiguity of the sign before the radical is determined from the following 
considerations. When /9 is a somewhat large positive quantity it is evident 
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that we ought to take the upper sign, as always approximately rf^ = o!'^, and 
then n} ^= 4l(i'^ + 4^8, thus n^ is larger than 4n'^. Hence, if we begin with n^ 
larger than 4w'^ and keep diminishing ^, continuity demands that the upper 
sign be maintained. And we may persist in diminishing /3 until, if 2e^ — e'' is 
a positive quantity, ^ ^= +T/6'^(2e^ — e'^) when n'^ =^ a'^ — \y'b''^{^^ — e''') and 
n^ ■= 4a'^ ; thus n^ is still larger than 4?(.'^. But if /3 is a somewhat large nega- 
tive quantity it is evident that the lower sign must be attributed to the radical ; 
and then approximately w'^ = a'^, n^ = 4«'^ + 4/3 ; that is, here n^ is less than 
4«'^. Then /3 may be numerically diminished until /3 = — i/b'' (2e" — e'^), when 
7t'^^ «"-^ + ^T/6'^(2e'^ — e'^) and «^ = 4a'^ ; and n^ is still less than 4n'l On 
the other hand if 2«^ — e'^ is a negative quantity, /9 in both cases can be nu- 
merically diminished until it vanishes, and then, in the first case, we have 
w'2 = a'2 + ^i/J'2 (e'2 _ 2e^\ ri" = 4a'2 + 4]/ 6"^ {e"" — M), and, in the second 
case, n'2 = a'^ _ ^y i'-^ (e"" ~ ^e") , ii' = 4«'2 — 4i/S'2(e'2 -~2?). 

From all this it is apparent that we cannot have 4?t'^ — w- ^ unless two 
conditions are satisfied. The first is /3 ^ 0, and the second either b' ^ or 
2^ — e'^ = 0. If the second condition is not fulfilled the limits towards which 

1 7' '2 7.' ' 

tend the coefficients ^ -= ^-7= and -i= y^„ , severally belonging to cos 

2nr — 4w^ n^ — {n — nf •' ° ^ 

2(p' in X and cos ((f — (f') in a;', can be written so as to cover all the cases 

, 1 6-2 , ee 

and 



- 4 1 / ± (2e2 — e'2) - ,/ + ^2e^ _ g'2) ' 

The ambiguous sign within the radical must be taken so as to render the quan- 
tity following positive ; that without the radical depends at once on the sign 
of b' and on whether /3 has been supposed to approach the limit from positive 
or negative values. These limits are independent of b', and it will be perceived 
that instead of being of the second order with respect to e and e' they are of 
the same order as the first terms of x and so'. The second supposition of the 
second condition viz. 2«^ — e'^ = renders these limits infinite ; but it is prob- 
able that this conclusion would be modified if terms of the fourth order were 
included in the equations determining n and n'. 

There is no need of giving the similar investigation for the assumption 
n' =- 2«. All the results of the latter are obtained from those of the former 
assumption by simply interchanging the accents. Thus the limits of the co- 
efficient of cos 2^ in x' and the coefficient of cos {f — f') in x are severally 

_i_ 1 ^ d -4- - ^-^^' 

~ 4 i/±(2e'2=^^) ^^ ~ |7'"±°(2F"— "?) ■ 
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In all this we have assumed that the quantity under the radical sign in 
the expressions for tn? and n'^ must not be negative. This is simply to insure 
stability of motion. If the constants of the problem are such as to make it 
negative, or if either v? and n'^ come out negative, the problem has still a real 
solution but the motion is unstable. 

The question of periodic solutions is intimately connected with the van- 
ishing of integrating divisors. If we suppose that in -\- i'n' = 0, where i and 
i' are integers prime to each other, by adopting an argument ([> such that 

-Vt = i' -n and — ^ = — i -^ , the series for x and x' reduce to the form 
at at dt at 

a; or a;' = Oq + a.^ cos ^ + /3i sin (p -\- a^ cos lip -\- ^^ sin 2^ -j- • • • • 

The variables then return to the same values after ^ has gone through a cir- 
cumference as also do their differential coefficients. If we conceive a system of 
rectangular coordinates for graphically exhibiting the values of x and x', the 
representative point P, in its motion, will describe a certain curve which, if the 
variables are confined within limits, will lie within a limited portion of the 
plane. This curve may have multiple points ; let us see what is the condition 
necessary and sufficient that this may have place. The motion of x and x' 
being a constant oscillation between minimum and maximum values, the curve 
must be either a spiral having no multiple points, in which case there would 
be a progressive diminution or augmentation of the minimum and maximum 
values of x and x', or there must be an infinite number of multiple points. The 
periodic solution is the case which lies between these two classes of solutions. 
We may conceive the angle of intersection at the multiple point to constantly 
diminish through variation of the arbitrary constants of the problem until it 
vanishes, when the periodic solution is reached. 

Let us suppose that the solution of the system of differential equations 



admitting the integral 



(Px ^ 3il (ZV ^ 3ii 

cie ~ dx ' de ~ 9x" 



1 dx^ -f- dx'^ f, , ri 



has been discovered and that it may be written 

x=.i^(t), x' = {^{t). 

Let t' deno'te another value of the independent variable t ; then, if the equa- 
tions 
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t and f being regarded as the unknowns, are satisfied by real values of t and tf, 
the curve under discussion will have multiple points ; but, if the roots of these 
equations are all imaginary, the curve will be a spiral. It appears that when 
X and a;' are contained within finite limits the roots t and t' are all real and 
infinite in number ; and thus that, in general, there are an infinite number of 
multiple points. 

Let us now suppose that the arbitrary constants are adjusted so that, in 
addition to the equations f, (f) = t^ {() . fj {(') = i^ (t) being satisfied by certain 
values of t and <', the equation 

W) fi'(0 

is satisfied by the same values. (Accents have been used to denote differen- 
tiation of the form of f.) Then it is evident that the last equation is the con- 
dition of the presence of a periodic solution ; and if we put t' — t = P, the 
equations 

f.(^ + P) = fi(0> i,{t^P) = i,{t), 

are satisfied for all values of t. 

A periodic solution must have at least one arbitrary constant less than in 
the general case. But, for the differential equations we are treating, there is 
usually a periodic solution in which there are only two arbitrary constants. 
This happens when both variables arrive simultaneously at their maximum or 
minimum values. Integrating the differential equations by the application of 
Maclaurin's Theorem, we should have, in general, 

«,==« + j^ + ^_ + ..., 

where a, a', b, h', are the arbitrary constants, and in the derivatives of ii it is 
understood that x and a/ are to be replaced severally by a and a' . In the case 
where both variables arrive simultaneously at a maximum or minimum value, 
and we count t from this epoch, not only J =^ 0, J' = 0, but also all the coefficients 
of the odd powers of < in a; and a?', as is easy to see from the successive differ- 
entiation of the derivatives of ii with respect to the time. We may therefore 
put 

a; = a -j- a^p + a/ + aji" +..., = f, {t) , 

a;' = a' + d^e + dj^ + a\e + =iS) ■ 
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In the previous equations for determining the times of the representative point 
P passing through multiple points, we may suppose <= — t, and the equations 
fi( — = fi(0) U.{.— = 'sW ^6 fulfilled for all values of t. Also the equa- 
tion showing the presence of a periodic solution is, in this case, an identity. 
This is explained by the circumstance that the point P having attained the 
limiting curve whose equation is i2 + C = returns on its previous path. 
Hence the condition necessary for the existence of a periodic solution in this 
case is that when the path is suflBciently prolonged backward it should again 
meet the limiting curve for real velocities. The motion of P will then be a 
swinging back and forth on this arc. This is tantamount to saying that the 
values of a and d must be so adjusted that a value t^ for t can be found which 
will make x and a;' have such values that they satisfy the equation U -\- (7=0. 
It may be of interest to know the earlier terms of the expansions of x and 
a;' in powers of ^. For the sake of brevity we put 



2^=r^r+ 



Then 



ySx^ 



'dQV 
9x'j 






9x1.2 ' dx 1.2.3.4 
'^ I \_Sx Sx" "^ 9x' 5xSx'\ \_3a? dx '^ 3xSx' Sai \\ &\^ ' " 

'^ I \_9ii 3x'^ 3x 5xSx'\ \jx^ 9x' '^ 3xSx' 3x\ ^ 6! "^ ' " " 

where in the partial derivatives it is supposed that x and x' are replaced sev- 
erally by a and a'. 

Evidently, for these series proceeding according to ascending powers of 
<^, we may substitute others proceeding according to cosines of multiples of an 
angle proportional to t, or x and a;' may have the form 

a; or a;' = a„ + a, cos nt -\- a^ cos 2n< + . . . . 

But here we are not entitled to assume that the coeflScients rapidly diminish. 
To bring about such a relation as w = w' it may be necessary for the arbitrary 
constants we have denoted by e and e' to be of the zero order of magnitude. 

If in the expansion of Q according to powers and products of x and x' we 
omit all the terms involving both variables the integration of the differential 
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equations is reduced to quadratures. Let us write in a more general way this 
expansion of Q 

+ a^* -\- tfgi a^x' ~\- a^^x'^ + a^^ xx'^ + (tot"'* 



It must be noticed that the coefficients of the terms involving both x and x' 
are of the order of the planetary masses while those which multiply the simple 
powers of x or x' are of the order of the central mass. If we put 

i2o = «2oar' + a^aff^ + a^^ + . . . 

the differential equations resulting from this expression of the potential func- 
tion are 

-^ = %a^x + Sajo*" + 4a^«« + . . . , 

-^ = 2ff^' + 3aos«" + 4ao4a?" + 

From these are immediately derived the integrals 

J^ ~ ^ + ^a)*" + «30-*'' + «««* + • • • , 

dvP' 

-^ = C" + a^'^ + aos«" + ^o^a? -f • • • , 



and thence the final integrals 

, , _ r dx 



. + .'=/ 



dx! 



The polynomials under the radical sign, on being equated to 0, will have each 
two real roots, which, on account of the smallness of x and a;', will be, sever- 
ally, in the neighborhood of the values 

a- = +A > »; = —-»- , ■'' = +A — — ' ^ — ^-vr 
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Tlie exact values of these can readily be found by trial ; let them be denoted 
severally, by a, b, a', b'. Then we can write 

, _ r dx 

^ + "^ ~ J ,/ [(a ^x){x-h){Ao+ A,x + A^-^ + ...)] ' 

. , _ r ^ 

^"^ ~J y/ [(a' - x') {x' - b') {A\ + A\x' + A'^'-' ...)]• 

The limiting values of x are a and b, and those of «', a' and b'. Within these 
limits the last factors under the radical sign are always positive and do not 
vary much, and thus x and x' are periodic functions of t, the periods being 
given by the formulae 

y [(a -«) {X - b) (^0 + A,x^A^ + .. .)] 
dx' 



T' = 2f 



^[(a' - x) {x' - b') {A', + A\^ + A'.^x'^ + •••)] ' 
By adopting, in place of x and x!, variables (}' and tp' such that 

2a; = (a + b) — (a — b) cos <l> , 2a;' = (a' + b') — (a' — b') cos tp' , 

the integrations are considerably simplified. The motion of (p and ^' is unlim- 
ited like that of t. Let 

[^„ -^ A,x-\- A^' + . . .]-J =z^[l + B,coa^ + 2B, cos 25^ + ...], 

[^'0 + -4'i«' + ^>'2 + . . .]-i = ^ [1 + ^'1 cos f + 2B, cos 2<p' + ...]. 

Then 

n{t + e)=(p + Bi sin ^ + ^2 sin 2^ + . . . , 



n 



'{t + c') = ^' + i?'i sin ^' + B'^am ^ij/ + . . . . 



The values of ip and ^' can be derived from these transcendental equations 
and have the forms 

(p = n(t-{: c) + E^ cos w (< + c) + ^2 cos 2 (»< + c) + . . . , 
^' = n'{t + c') + E{ cos w' (!! + c') + jE^' cos 2 (?i'< -[- c') + . . • • 
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From these we can pass to the values of cos ^ and cos (p', and thence to those 
of X and x', which have the form 

X ^ i^ G^ -\- G^ cos f -\- G^ cos 2^ + . . . , 

a;' =: ^ G^ + G^ cos f' + G^ cos 2^' + . . . , 

^ and ^' being put for nt -\- c and n't -\- c'. 

The reversion of the periodic series can be avoided by the use of definite 
integrals, and we can pass directly to the values of the Gi and G^. For we 
have 



2 {* 
Gt^= - \ X cos i<p d<p , 



with a similar formula for Gi'. By adopting ^ as the independent variable 
this becomes 



2 /* 
6^4 = - I [(a + b) — (a — b) cos </>] cos [i(p + iB^ sin <p + e^j ^^^ ^'r' + • • -1 



[1 + ^1 cos ^ + 2B2 cos 25^ + . . .] <:?^ . 

In the particular case of i ^0, this reduces to 

6^0 = a + b - i (a - b) i?i . 

Proposing now to take into account the omitted terms of iJ, which we will 
denote as i2], let the values of x and x', just determined, be denoted as x„ and 
x„', their corrections as dx and dx', the differential equations for the latter will 
be 

-^ = [2«a, + Ba^{«! + %) + 4a«(«^ + a;a;„ + ajj^) + . . .] ^a; + -^ , 

^ = [2«„, + 3«^ (*' + x',) + 4a<^ (a;'^ + x' x', + ^V') + ...-\dx' + ^. 

These can be integrated by successive approximations, introducing into the 
right members the values of x and x' from the previous approximation, begin- 
ning with X = x^, 3! = x\. In order to prevent t from getting outside of the 
functional sign cos, we adopt the following device : Let n and n' be two in- 
determinate constants, whose values come out only at the end of the approxi- 
mation, but which are equivalent to the rigorous rates of motion of the 
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arguments <p and if'. We add n^dx to each member of the first equation and 
n'^dx' to each member of the second. Then, putting 

9Q 
Q=\n^ + %%^ + 3a,o(« + *o) + 4a«(«^ + ««„ + x^) + . . .] ^a; + -^, 

and ^ for the similar quantity in the second equation, we have 

^ + n^bx =q, ^ + «'W = ^ . 

In any stage of the approximation we shall have 

Q^ I Ki^i. cos {i(p + i'(p') , Q' = I K\t, cos {if + i'(p'} , 
from which 

dx ^ 2'. [w'* — (en + i'n'f]^^ iLif. cos (2^ + i'<p') , 

Then the equations for determining n and n' are 

-^1,0 = > -^'0,1 "^ ^ • 

The terms of i<? in ar^a?' and xx'^ can be removed from it by a linear trans- 
formation, but at the expense of reintroducing a term in xx', as also a term 
involving the product of the velocities into the expression for the living force. 
For substituting 

« = y + V > «' = y' + i'y , 

if I and I' are so determined that 

da J + a,, a + m) + a J (2 + W) + ZaJ' = , 
da J' + a J (2 + W) + a,2 (1 + 2«') + ZaJ' = , 

the terms in U involving ify' and yy''^ vanish. And if we put 

-"20 = ^20 + ^n'' ~r ''W ' 

5„ = 'ia.J + «!, (1 + IV) + 2a J' , 

•O02 = '■^20^ ~l" ^11^ "l~ '''02 ) 

-030 = «30 + ^21'' + <''l2^'^ + <^0S*' ) 

-Oo3 :="- <<3o^ "h <*2/ ~r "^12^ + '''os ) 
we have 

a = B^f + A,3/3/' + B^' + i^soS/" + B^y'' + . . . , 
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and the expression of the living force becomes 

and the differential equations of the problem are 

(1 + Z'^) ^ + (Z + I') ^ = IB^y + ZB^f + B,,y' , 

(1 + ^')^ + (^ + ^ = 2^0^ + ZB^y'^ + B,,y . 

The elimination of I or I' between the equations which determine them 
gives rise to an equation of the fifth degree, which consequently, always has a 
real root. As ajj and a,^ are of the order of the masses of the planets, ap- 
proximate values for I and I' are I = — ~^ , ^' = — o^ • 

In integrating the preceding equations one may follow either of two 
courses ; first, we may neglect the terms in y' and y'^, which will afterwards 
be regarded as perturbative, and then we shall have linear equations with con- 
stant coefiScients easily integrable ; or second, we may regard the last terms of 
both members of each equation as perturbative, and the equations restricted 
to the remaining terms are integrable by quadratures, as has been shown. 



